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SPECIFICATION 
Method and Equipment for Measuring Feature Points of Wave Sig nal 
FIELD OF THE INVENTION 

The present invention relates to a method and equipment for detecting with 
5 excellent accuracy the feature points of wave signals with irregular feature point 
values on the waveform or irregular distance between feature points, and more 
specifically, relates to a method and equipment for measuring feature points of wave 
signals which can be ideally applied to the measurement of such objects as the 
number of tree rings in a piece of wood or the width of the tree rings. 

10 BACKGROUND ART 

In the field of dendrochronology, it is possible to establish in annual units 
the year in which each of the tree rings in a given piece of wood was formed, by 
cross-referencing against a database of standard tree ring width fluctuations. This 
serves as the basis for tree ring dating. This database, brought about by the 

15 intensive efforts of the National Research Institute for Cultural Properties, Nara, now 
enables researchers in Japan to go back to 912 B. C. for hinoki cypress, and to 1313 
B. C. for sugi cedar. Incidentally, in Germany, a nation that is at the forefront of 
dendrochronology, standard databases have been created with a span of 
approximately 10,000 years. The field of dendrochronology deals primarily with 

20 the following matters: 

(a) Estimation of the year of felling of a piece of wood 

(b) Estimation of the year of creation and course of repair of wooden 
cultural properties (architecture, Buddhist carvings, works of art and handicrafts, 
etc.); authentication; etc. 

25 (c) Study of climatic changes over long periods of time in the past; 

study of global warming; etc. 

The ultimate in detection performance is required of the time series data of 
each tree ring width that is used in dendrochronology: both erroneous 
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detection(erroneously recognizing something that is not a tree ring as a tree ring) and 

non-detection(failing to recognize a tree ring) must be zero. 

For this reason, the measurement of tree ring width has been carried out by 

the human eye, via specialized systems that use a measurement microscope. Such 
5 work has required a high level of skill and enormous amounts of time (approximately 

1 hour for a specimen with 300 or so tree rings). The large scale of the system 

setup was another problem inherent in this method. 

In an effort to automate the measuring work, methods have been considered 

wherein each tree ring width is measured using a personal computer to analyze tree 
10 ring images acquired using such image acquisition equipment as digital cameras and 

scanners. Several endeavors have been made along these lines to date, but the 

current situation is one in which problems such as the aforementioned detection 

performance requirements and large scale of the system, as well as price 

considerations, have kept such methods from becoming widespread as a means for 
15 research. 

In particular, the problem with Japanese cypress (Hinoki) has been that, 
despite its wide use in cultural properties and hence its importance as a 
dendrochronological species, restrictions such as the narrowness of tree ring width 
and the indistinctness of tree rings compared to Japanese cedar (Sugi) cedar have 
20 made the practical application of automated measurement extremely difficult. The 
following are the major publicly known technologies that are similar to the present 
invention: 

(1) "Win DENDRO," Regent Instruments, Canada 1988 
(see http://www.regent.qc. ca.products/dendro/DENDRO.html) 
25 Designed by Dr. Rejean Gagnon and Dr. Hubert Morin of Quebec University 

and commercialized by Regent Instruments, this software was developed for 
dendrochronological research. This software is able to conduct tree ring 
measurement and wood tissue analysis on the basis of information on light intensity 
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differences in the tree ring image. While the details of the algorithms of this 
software are unknown, as far as can be surmised from the wording in the company's 
catalogue, said software does not appear to use wavelet processing or technology to 
integrate information from multiple measuring lines. 
5 (2) "Gazo shori shisutemu wo mochiita nenrin haba keisoku (Measurement 

of tree ring width using an image processing system)," Noda, Masato 1990: 
Presentation at the Tree Ring Society 

This presentation concludes that, while it is possible to measure the tree 
rings of Japanese cedar (Sugi), it is impossible to measure those of Japanese cypress 
10 (Hinoki). Measurement methods relating to this presentation do not use wavelet 
processing or technology to integrate information from multiple measuring lines. 

(3) Japanese Unexamined Patent Publication(Kokai) No. HI 1 - 232427 

There is a description of the use of light intensity information in the image to 
measure the number of tree rings; however, said technology is already publicly 
15 known due to (2) above. Neither wavelet processing nor technology to integrate 
information from multiple measuring lines is used in any way in the publicly known 
technology listed in this publication. 

By acquiring pixel information from the tree ring image along a 
measurement line, it is possible to obtain waveform signals of information on light 
20 intensity changes and/or waveform signals of information on density changes. The 
maximum point of the density waveform (or in the case of the intensity waveform, 
the minimum point) corresponds to the darkest portion(the highest density late wood 
portion) of each tree ring layer. Therefore, by recognizing the maximum point of 
the density signal waveform or the minimum point of the light intensity waveform, it 
25 is possible to recognize each tree ring layer. 

On the other hand, further treatment of the density waveform signal by 
differential processing makes the dark to light transition point (the minimum point of 
the differential waveform) correspond to the end point (late wood end) of each tree 
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ring layer. Therefore, by measuring the distance between minimum points of the 
differential waveform signal on the measurement line, it is possible to measure tree 
ring width with greater accuracy. 

The obtainment of waveform signals of information on light intensity 
5 changes and/or waveform signals of information on density changes by acquiring 
pixel information from the tree ring image along a measurement line is a publicly 
known matter due to the publicly known literature described above. 

However, the tree ring widths of wood specimens are generally irregular, 
and it is not unusual to encounter up to 100-fold differences between the maximum 

10 tree ring width and minimum tree ring width. Therefore, when the detection 
accuracy for the feature points (the peak points, which are the maximum points, or 
the trough points, which are the minimum points) for the waveform signal obtained 
from the tree ring image is set to the level of detecting the small distances between 
feature points, the analysis is prone to picking up noise unrelated to tree rings in 

15 portions where the distances between feature points are large. On the other hand, 
when the detection accuracy for the feature points (the peak points or the trough 
points) for the waveform signal obtained from the tree ring image is set to the level 
of detecting the large distances between feature points, the analysis may fail to detect 
feature points. Therefore, it is difficult to measure with accuracy the number and 

20 width of tree rings in specimens with large differences between the minimum tree 
ring width and maximum tree ring width. 

Furthermore, because the density level is not uniform in the tree ring image, 
it is often the case that the light intensity waveform signal, density waveform signal, 
differential waveform signal, etc. obtained from the tree ring image all have 

25 undulating features over the entire interval to be measured. For this reason, 
attempting to use a fixed threshold value to detect feature points (the peak points, 
which are the maximum points, or the trough points, which are the minimum points) 
in the waveform signal can result in a failure to detect feature points, leading to an 
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inability to measure the number of tree rings or tree ring width with accuracy. 

Therefore, a measurement method and equipment that allows for the speedy 
and highly precise acquisition of time series data on each tree ring width, which is 
the most basic data in the study of dendrochronology, is desired. 
5 DISCLOSURE OF THE INVENTION 

In order to provide a method and equipment for measuring feature points of 
a waveform signal which can comply with the aforementioned desire, the present 
invention aims to provide a measurement method and equipment capable of speedy 
and high precision detection of waveform feature points even when a waveform 

10 signal has irregular feature point values or irregular distances between feature points. 

In order to resolve the above-mentioned problems, the present invention, as 
the first invention, provides a method for measuring feature points of a waveform 
signal having irregular feature point values or irregular distances between the feature 
points, the method comprising the steps of : performing wavelet conversion of a 

15 waveform signal within a predetermined interval by using a predetermined mother 
wavelet and multiple scale levels; calculating squared mean for interval for each 
interval width corresponding to said scale levels in relation to a wavelet conversion 
signal for each scale level generated by the said wavelet conversion; defining a scale 
level at a point where the calculated value of the said squared mean for interval 

20 becomes maximum at an arbitrary point within the predetermined interval, as a 
dominant level for that point; and detecting points at which the said waveform signal 
reaches maximum value or minimum value for each interval width corresponding to 
the dominant level, as the feature points of the waveform signal. 

Furthermore, as the second invention, in the measurement method having the 

25 constitution of the above-mentioned first invention, the present invention provides a 
method for measuring feature points of a waveform signal, wherein the 
aforementioned wavelet conversion uses the following formula (6), that is, 
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dj(x) = b J [j(b J (x - k))f(x)dx 



(6) 



where f(x) is the waveform signal, \|/(x) is the mother wavelet, W is a scaling 
parameter, b is a constant (b > 1), j is a scale level comprised of zero or a negative 
whole number, and k is a translating parameter. 

Furthermore, as the third invention, in the measurement method having the 
constitution of the above-mentioned second invention, the present invention provides 
a method for measuring feature points of a waveform signal characterized in that the 
aforementioned mother wavelet uses a French hat wavelet transform which is defined 
by the following formula (7), that is, 



0(x) = 



1 - 1 < X < 1 
-0.5 -3<x<-l, or \<x<3 -(7) 
0 x<-3 } or 3 < jc 



<K*) = — — e ~ X =0-2*V ...(8) 



Furthermore, as the fourth invention, in the measurement method having the 
constitution of the above-mentioned second invention, the present invention provides 
a method for measuring feature points of a waveform signal characterized in that the 
aforementioned mother wavelet is a Mexican hat wavelet transform which is defined 
15 by the following formula (8), that is, 

2dx" 

Furthermore, as the fifth invention, in the measurement method having the 
constitution of the above-mentioned second invention, the present invention provides 
a method for measuring feature points of a waveform signal characterized in that the 
20 aforementioned calculation of the squared mean for interval uses the following 
formula (9), that is, 

JX-pj 



s,(x) = 2- , p; , r''h<*)f <ft 



(9) 
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where j is the scale level used in the formula (6), k is the translating parameter, and 
Pj is a constant that is set according to scale level j so that the constant pj becomes 
larger as the scale level j becomes lower. 

Furthermore, as the sixth invention, in the measurement method having the 
5 constitution of the above-mentioned fifth invention, the present invention provides a 
method for measuring feature points of a waveform signal characterized in that pj in 
the aforementioned formula (9) for the calculation of the squared mean for interval is 
defined by the following formula (10), that is, 

Pj=b~ J a ...(io) 

10 where a is a constant determined by the support of the mother wavelet \|/(x) used in 
the formula (6), b is the constant used in the formula (6), and j is the scale level used 
in the formula (6). 

Furthermore, as the seventh invention, in the measurement method having 
the constitution of the above-mentioned second invention or fifth invention, the 

15 present invention provides a method for measuring feature points of a waveform 
signal characterized in that the value of b in the aforementioned formula (6) is 2. 

Furthermore, as the eighth invention, in the measurement method having the 
constitution of the above-mentioned first invention, the present invention provides a 
method for measuring feature points of a waveform signal characterized in that the 

20 aforementioned waveform signal is a pixel light intensity or density information 
signal acquired from a target image, such as wood specimen tree ring image or the 
like, along a measurement line configured on the image. 

Furthermore, as the ninth invention, in the measurement method having the 
constitution of the above-mentioned first invention, the present invention provides a 

25 method for measuring feature points of a waveform signal characterized in that the 
aforementioned waveform signal is a pixel light intensity or density information 
signal acquired from a target image, such as wood specimen tree ring image or the 
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like, along a measurement line configured on the target image that is further 

subjected to differential processing. 

Furthermore, as the tenth invention, in the measurement method having the 

constitution of the above-mentioned ninth invention, the present invention provides a 
5 method for measuring feature points of a waveform signal characterized in that the 

aforementioned differential processing is a calculus of finite differences between 

multiple pixels separated by an interval of several pixels. 

Furthermore, as the eleventh invention, in the measurement method having 

the constitution of the above-mentioned eighth invention, the present invention 
10 provides a method for measuring feature points of a waveform signal characterized in 

that when the aforementioned waveform signal is a density information signal, said 

density information signal is f(x), the aforementioned dominant level is jd, the 

constant corresponding to said dominant level is qjd, and an arbitrary point on the 

aforementioned measurement line is x, then when the value of f(x) is equivalent to 
15 the maximum value f ma x(x) of f(x) of the interval [x - qj d , x + qjd], the point x is 

determined as the feature point which indicates the maximum density point within 

the tree ring layer. 

Furthermore, as the twelfth invention, in the measurement method having the 
constitution of the above-mentioned ninth invention, the present invention provides a 

20 method for measuring feature points of a waveform signal characterized in that when 
the aforementioned waveform signal is a differential signal obtained by differential 
processing of a density information signal, the said differential signal is f(x), the 
aforementioned dominant level is jd, the constant corresponding to said dominant 
level is qjd, and an arbitrary point on the aforementioned measurement line is x, then 

25 when the value of f(x) is equivalent to the minimum value f m j n (x) of f(x) of the 
interval [x-qj d , x + qjd], the point x is determined as the feature point which 
indicates the late wood end within the tree ring layer. 

Furthermore, as the thirteenth invention, in the measurement method having 
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the constitution of the above-mentioned eighth invention or ninth invention, the 
present invention provides a method for measuring feature points of a waveform 
signal characterized in that the aforementioned measurement line is comprised of a 
main measurement line and multiple subordinate measurement lines which are 
5 equidistant parallel lines on either side of said main measurement line, and when 
waveform signal feature points are detected at a point that is the same distance from 
the starting end on said main measurement line and subordinate measurement lines, 
then those feature points are determined to be a feature point on the main 
measurement line provided that one of the conditions is that the number of said 

10 feature points comprises at least a majority in relation to the number of measurement 
lines including the main measurement line and subordinate measurement lines. 

Furthermore, as the fourteenth invention, in the measurement method having 
the constitution of the above-mentioned thirteenth invention, the present invention 
provides a method for measuring feature points of a waveform signal characterized in 

15 that two subordinate measurement lines are configured respectively at both sides of 
the aforementioned main measurement line, and when feature points of the waveform 
signal are detected at a point that is roughly the same distance from the starting end 
on the said main measurement line and subordinate measurement lines, and when 
feature points are found on the main measurement line and on at least one of the two 

20 subordinate measurement lines that are positioned adjacent to said main 
measurement line, and when feature points are found on the two subordinate 
measurement lines that are positioned adjacent to the main measurement line and on 
at least one of the other subordinate measurement lines, then those feature points are 
determined to be a feature point on the main measurement line. 

25 Furthermore, as the fifteenth invention, in the measurement method having 

the constitution of the above-mentioned eighth invention or ninth invention, the 
present invention provides a method for measuring feature points of a waveform 
signal characterized in that a smoothing process using peripheral pixel information is 
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performed on the pixel light intensity or density information acquired from the target 
image along the measurement lines configured on the target image. 

Furthermore, as the sixteenth invention, the present invention provides an 
equipment for measuring feature points of a waveform signal having irregular feature 
5 point values or irregular distances between feature points, characterized by 
comprising: a wavelet conversion means for performing wavelet conversion of a 
waveform signal within the predetermined interval by using a predetermined mother 
wavelet and multiple scale levels; a squared mean calculation means for calculating 
squared mean for interval for each interval width corresponding to said scale levels 

10 in relation to the wavelet conversion signal for each scale level generated by the 
wavelet conversion means; a dominant level decision means for defining a scale level 
at which the calculated value of the aforementioned squared mean for interval 
becomes maximum at an arbitrary point within the aforementioned predetermined 
interval, as the dominant level for that point; and feature point detecting means for 

15 detecting points at which the aforementioned waveform signal reaches maximum 
value or minimum value for each interval width corresponding to the dominant level, 
as the feature points of the waveform signal. 

Furthermore, as the seventeenth invention, in the measurement equipment 
having the constitution of the above-mentioned fifteenth invention, the present 

20 invention provides an equipment for measuring feature points of a waveform signal 
characterized by further comprising distance calculating means for calculating 
distances between the feature points on the basis of the detected feature points of the 
waveform signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 FIG. 1 is a flow chart showing the first half of the process of the method for 

measuring feature points of a waveform signal according to one embodiment of the 
present invention applied to tree ring measurement. 

FIG. 2 is a flow chart showing the second half of the process that follows the 



flow chart in FIG. 1. 

FIG. 3 is an explanatory diagram showing typically an example of the 
configuration of the main measurement line and the subordinate measurement lines. 

FIG. 4A is an explanatory waveform diagram of a RGB signal showing the 
5 manner in which the RGB signal acquired along the measurement line is converted 
into a density signal, FIG. 4B is an explanatory waveform diagram of a density 
signal generated from the RGB signal, and FIG. 4C is an explanatory diagram of the 
RGB -density conversion curve used in the conversion of the RGB signal into the 
density signal. 

10 FIG. 5A is an explanatory waveform diagram showing the changes in the 

density signal, and FIG. 5B is an explanatory waveform diagram showing the density 
signal after smoothing. 

FIG. 6A is an explanatory waveform diagram of the density signal, FIG. 6B 
is an explanatory waveform diagram of the differential signal generated by 
15 differential processing of the density signal, and FIG. 6C is an explanatory diagram 
for explaining the method of differential processing. 

FIG. 7 is an explanatory waveform diagram showing an example of an 
extension dummy signal added to both sides of the density signal. 

FIG. 8A is an explanatory diagram showing a typical example of the French 
20 hat type mother wavelet used in wavelet conversion, and FIG. 8B is an explanatory 
diagram showing a typical example of a Mexican hat type mother wavelet. 

FIG. 9 is an explanatory diagram showing the manner in which the French 
hat type mother wavelet waveform changes in stages according to the level value. 

FIG. 10 is an explanatory diagram showing the signal waveform that is 
25 obtained by wavelet conversion of the density signal using three graded level values. 

FIG. 11 is an explanatory diagram showing the results of the processing of 
squared mean for interval according to the level value applied to the signal waveform 
obtained by the wavelet conversion shown in FIG. 10. 
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FIG. 12 is an explanatory diagram showing the results of the decision 

i 

operation of the dominant level from the results of the processing of the squared 
mean for interval shown in FIG. 1 1 . 

FIG. 13A is an explanatory diagram of a dominant level showing 
5 schematically the process of searching and deciding the maximum value within the 
interval on the density signal waveform and the minimum value within the interval of 
the differential signal waveform, using the results of the dominant level determinant 
operation shown in FIG. 12, FIG. 13B is an explanatory waveform diagram of a 
density signal, and FIG. 13C is an explanatory waveform diagram of a differential 
10 signal. 

FIG. 14A is an explanatory diagram showing the results of the detection of 
the maximum density point (feature point) from the density signal, and FIG. 14B is 
an explanatory diagram showing the results of the detection of the late wood end 
point (feature point) from the differential signal. 
15 FIG. 15 is an explanatory diagram showing schematically the process of 

establishing a detection point (feature point) by cross-referencing of the main 
measurement line and subordinate measurement lines. 

FIG. 16 is an explanatory diagram showing schematically the process of 
establishing a detection point (feature point) by cross-referencing of the subordinate 
20 measurement lines following the process in FIG. 15. 

FIG. 17 is an explanatory diagram showing schematically the final result of 
the process of establishing a detection point (feature point) by cross-referencing 
following the process in FIG. 16. 

BEST MODE FOR CARRYING OUT THE INVENTION 
25 Preferred embodiments of the present invention will now be described with 

reference to the drawings. 

FIG. 1 and FIG. 2 are flow charts showing one embodiment of the method 
for measuring feature points of a waveform signal according to the present invention 



13 

applied to a tree ring measurement. Referring to FIG. 1 and FIG. 2, the process of 

i 

tree ring measurement is comprised of the following steps: 

1. Acquisition of Tree Ring Image (Step SI) 

2. Designation of Measurement Site (Step S2) 
5 3. Acquisition of Pixel Information(Step S3) 

4. Density Conversion(Steps S4 and S5) 

5. Smoothing (Steps S6 and S7) 

6. Acquisition of Peak Signal or Edge Signal (Steps S8 through SI 2) 

7. Designation of Mother Wavelet (Step S13) 
10 8. Wavelet Conversion(Step SI 4) 

9. Calculation of Squared Mean for Interval (Steps SI 5 and SI 6) 

10. Determination of Dominant Level (Step SI 7) 

1 1 . Determination of Maximum Value within Interval (in the case of edge 
signal, the minimum value within the interval) (Steps SI 8 through S20) 

15 12. Cross-referencing of Main Measurement Line and Subordinate 

Measurement Lines (Step S21) 

13. Cross-referencing between Subordinate Measurement Lines (Step S22) 

14. Determination of Tree Ring Location(Step S23) 

15. Measurement of Tree Ring Width and Output of Results (Steps S24 and 

20 S25) 

Furthermore, the equipment for measuring feature points of a waveform 
signal according to one embodiment of the present invention is a program which is 
capable of executing the above process from step S2 through step S25, and becomes 
capable of execution by installing into a computer device such as a personal 
25 computer into which an operating system has been installed. 

The program in this embodiment comprises: means for designating a 
measurement site, means for acquiring pixel information, means for performing 
density conversion, means for smoothing, means for acquiring peak signal or edge 
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signal, means for designating mother wavelet, means for performing wavelet 
conversion, means for calculating squared mean for interval, means for deciding 
dominant level, means for deciding maximum value within the interval (in the case 
of edge signal, the minimum value within the interval), means for cross-referencing 
5 of main measurement line and subordinate measurement lines, means for cross- 
referencing between subordinate measurement lines, means for deciding tree ring 
location, means for measuring tree ring width, and means for outputting results. 
Next, each of the above processes will be described in detail. 
1. Acquisition of Tree Ring Image (Step SI) 

10 The tree ring image is acquired using imaging equipment such as a digital 

camera, scanner, or the like. For example, specimens which are small and portable 
can be imaged by scanner, while large specimens can be imaged by high resolution 
digital camera. Here it is most important to pay attention to obtaining the required 
image resolution. Hinoki cypress specimens may present tree ring widths that are 

15 as narrow as approximately 0.1 mm per layer at places. In order to distinguish 
these portions, according to the Nyquist Sampling Theorem the size per pixel should 
be set at or under 0.05 mm, which is half of that width. 

Furthermore, in consideration of safety and measurement margin of error, it 
is preferable to conduct image acquisition at a resolution of 1,200 dpi, which is 

20 equivalent to the approximately 0.02 mm size of each pixel. In order to maintain 
this resolution over the entire measurement area of a large-size specimen, it is 
preferable to use the method of limiting the imaging area of each frame, and 
connecting the frames together. 

The tree ring images handled in the present invention can be, in addition to 

25 images obtained directly from the specimen surface using digital image acquisition 
equipment as described above, images of film photographs (photographs made using 
regular film), etc. that have been digitally converted. In addition, tree ring images 
inside the specimen can also be used, such as x-ray photographs, x-ray CT images (x- 
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ray tomograms), and MRI images. 

2. Designation of Measurement Site (Step S2) 

The measurement of tree ring width is conducted, as shown in FIG. 3, along 
at least one measurement line 1. On the cross section(transverse plane) disc, the 
5 measurement line is usually set from the pith (center left in FIG. 3) to the bark 
(center right in FIG. 3). On the radial section(longitudinal plane), the measurement 
line is set from the pith to the bark in a manner that intersects orthogonally with the 
grain. In the present invention, preferably, in order to improve the recognition 
performance, multiple subordinate measurement lines 2, 3, 4, and 5 are set on both 

10 sides of the single main measurement line 1, parallel to the main measurement line 
and spaced appropriately (see FIG. 3). 

The appropriate number of subordinate measurement lines, according to our 
studies, has been confirmed to be two for the cross section and four for the radial 
section. However, given individual differences among specimens, this number is 

15 not necessarily appropriate for all cases. On the actual program, it is preferable to 
designate the spaces and numbers of the subordinate measurement lines, as well as 
the starting point la and termination point lb of the main measurement line 1 on the 
tree ring image, to allow designation of the measurement area. 

3. Acquisition of Pixel Information(Step S3) 

20 As described above, the pixel information for each image (in the case of 

ordinary personal computers, the BGR 3 channel 8-bit digital signal) is acquired 
along the main measurement line 1 and the subordinate measurement lines 2 to 5. 
In relation to image acquisition, it is preferable to simultaneously acquire 
information on measurement line peripheral pixels for smoothing (the specific 

25 number of peripheral pixels is determined according to the level of smoothing) as 
described below, in addition to the pixel information on the measurement lines. 

It should be noted that, since the tree ring recognition by the methods of the 
present invention does not use color information, it is possible to use averaged 
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signals or mixed signals of BGR signals, or signals based on the light intensity of 

i 

signals of certain channels only (for instance, G signal only). Ordinarily, it suffices 
to use a signal that is an averaged mixture of B:G:R = 1:1:1. 
4. Density Conversion(Step S5) 
5 The pixel information acquired by the above method has the characteristic of 

being brighter as the value becomes greater, because the information is a result of the 
light intensity of each pixel. Furthermore, the RGB signal on the personal 
computer system has usually been processed via a nonlinear transformation called 
gamma correction in order to correct for the characteristics of the CRT monitor. On 

10 the other hand, since the light and dark contrast of tree rings is a result of the density 
of cells in the tree rings, it is usual practice in dendrochronology to describe light 
and dark contrast on the basis of density or image density which is related to density. 

The description of density becomes darker as the value becomes greater. In 
the embodiment according to the present invention, this usual practice of 

15 dendrochronology is followed, and the red and green averaged signal (RGB averaged 
signal) is converted to density signal (see FIGS. 4A, 4B, and 4C). For the 
conversion method, the conversion formula derived from the known density portion 
of the image (usually step tablet chart) can be used. One example of the conversion 
formula is as follows: 

D = -0.79 x\n(d) + 0.49 —(11) 

when : D is a dens i ty s igna I va lue f and 
d is a RGB averaged s igna I va lue 

20 

In addition, a linear conversion such as the following may be used in place of a non- 
linear conversion. 
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D = -fid + c -(12) 
when : J3 is a pos it ive constant va lue , and 
c is a constant va lue 

It should be noted that when conducting such image processing as binary 
coding, or when using conventional methods which necessitated the setting of a 
threshold value for tree ring recognition, the manner in which the tree ring signal 
5 waveform and amplitude characteristics were described was important. 

However, as described below, since in the present invention tree ring 
recognition can be conducted without regard to the tree ring signal waveform value 
itself, there is no great difference in tree ring recognition performance whether the 
signal used is RGB, light intensity, or density. Therefore, in order to shorten the 
10 processing step or time, this process can be omitted. However, it should be kept in 
mind that omitting this process will cause the tree ring signal waveform to invert. 

5. Smoothing (Steps S6 & S7) 

Smoothing is conducted using information of peripheral pixels on each of 
the pixel values obtained via the above process on the main measurement line 1 and 

15 subordinate measurement lines 2 to 5. By conducting smoothing, it is possible to 
reduce erroneous recognition attributable to noise resulting from the specimen itself 
or from the measuring equipment (see FIGS. 5A and 5B). 

As the method of smoothing, the moving means or moving median may be 
used. As the number of peripheral pixels, for example in the case of 1,200 dpi 

20 resolution, in our experiment the following were appropriate: cross section, approx. 5 
pixels in the direction orthogonal to the measurement line and approx. 3 pixels in the 
direction of the measurement line; and radial section, approx. 15 pixels in the 
direction orthogonal to the measurement line and approx. 5 pixels in the direction of 
the measurement line. However, these pixel numbers are not the only possibilities. 

25 It should be noted that when smoothing is not to be conducted, the smoothing range 
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is to be set as 1 pixel each in the direction orthogonal to the measurement line and 
the direction of the measurement line. 

6. Acquisition of Peak Signal or Edge Signal (Steps S8 through SI 2) 
In the signal acquired through the above process, the maximum point of the 
5 waveform represents the densest portion of each tree ring (late wood maximum 
density portion). This density information signal shall hereinafter be referred to as 
the peak signal. By recognizing the peak point (maximum point of waveform), 
which is the feature point of the peak signal, it is possible to distinguish each tree 
ring layer. 

10 On the other hand, the trough point (the minimum point of the waveform), 

which is the feature point of the differential signal (actually the signal obtained by 
calculus of finite differences) of the peak signal, corresponds to the point where the 
tree ring image shifts from the dark portion(high density portion) to the light 
portion(low density portion), and corresponds to the termination point of each tree 

15 ring layer (late wood end). This differential signal shall hereinafter be referred to 
as the edge signal. 

The ordinary differential signal is obtained by calculus of finite differences 
of adjacent pixels, but adjacent calculus of finite differences is prone to being 
affected by noise. Therefore, in the preferred embodiment of the present invention, 

20 efforts are made to reduce noise by conducting calculus of finite differences on 
multiple pixels spaced several pixels apart (see FIG. 6C). 

When the peak signal (the waveform shown in FIG. 6A) and the edge signal 
(the waveform shown in FIG. 6B) are compared, the peak signal has better detection 
performance for the identification of each tree ring, due to the effect of waveform 

25 noise, etc., but from the perspective of measuring tree ring width with accuracy, the 
edge signal, which corresponds to the late wood end, has higher precision. 
Therefore, in the preferred embodiment of the present invention, these characteristics 
are maximized by allowing the user to select from the following on a menu: 
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recognition by peak signal, recognition by edge signal, and recognition using a 
combination of both. 

In the wavelet conversion(convolution integration operation) described in 
the following section "8. Wavelet conversion" operations are conducted on pixels of 
5 a single point including peripheral pixels. Therefore, it is desirable for a dummy 
signal (see the dotted line portions on both ends of the waveform in FIG. 7) to be 
created on both ends of the measurement line (starting point la and terminating point 
lb in FIG. 3) by folding the signal back on itself, to prevent the operation from 
becoming inoperable due to lack of pixels. The length required for the dummy 

10 signal is dependent on the mother wavelet support designated below in "7. 
Designation of mother wavelet" and the level number designated in "8. Wavelet 
conversion", but since the longest length of end processing is at the lowest level j, it 
suffices to extend the dummy signal accordingly. 

7. Designation of Mother Wavelet (Step SI 3) 

15 The function \|/(x) that can be used as the mother wavelet is a function that 

meets the following two conditions: 

(a) The end portion of the function is 0 or converges on 0 

(b) The sum total of all intervals (integral value) is 0. 

In addition, it is desirable for the following condition to be met if possible, 
20 although it is not necessarily required: 

(c) The support (the interval at which the function value is not 0) is 

compact. 

As concrete example of a mother wavelet, the following French hat wavelet 
transform (formula (13)) and Mexican hat wavelet transform (formula (14)) can be 
25 given(see FIGS. 8A and 8B). 
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1 -1 <X < 1 

4)(x)= -1-0.5 — 3 < jc < -1, or \<x<3 
0 x<-3, or 3<x 



(13) 



2 <33C 

In the preferred embodiment of the present invention, the French hat wavelet 
transform of the formula (13) is used as default, because it is easy to generate and its 
5 support is compact. Of course the Mexican hat wavelet transform and other mother 
wavelets can be used if designated. 

8. Wavelet Conversion(Step SI 4) 

The wavelet conversion for the one-dimensional image signal f(x) is defined 
by the following formula (15). 

dj(k) = 2 J f 4>(2 J (x - k))f(x)dx ...(is) 

10 J-oo 

In the formula (14), f(x) is the peak signal or edge signal on the 
measurement line. In addition, as described in the above "6. Acquisition of Peak 
Signal or Edge Signal", both ends of the signal have been appropriately treated by 
folding back on the signal itself. v|/(x) is the mother wavelet described in "7. 

15 Designation of Mother Wavelet", 2 J is the scaling parameter, and the level (scale 
level) j is 0 or a negative whole number. By this scaling parameter, it is possible to 
change the width and height of the mother wavelet according to the power law of 2 
while retaining the area unchanged (see FIG. 9). k indicates the translating 
parameter. It should be noted that while in the formula (15) the scaling parameter 

20 is 2\ this is not the only possibility. It is also possible to have b* (where b > 1) and 
illustrate as in the following formula (16). 
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dj(x) = b j [j(b J (x - k))f(x)dx ...(16) 

In the formula (16), when 1 < b < 2 5 the spaces between the steps of the 
scaling parameter become finer than with the formula (14), and when 2 < b, the 
spaces between the steps of the scaling parameter become larger than in the formula 
5 (15). 

As has been verified by our experiment, a good detection performance is 
usually obtained with the power of 2 step represented in the formula (15). 

What is signified in the formula (15) and formula (16) is the convolution 
integration operation of and the mother wavelet \|/(x) which has had its size 

10 (stipulated by the scaling parameter) and position(stipulated by the translating 
parameter) changed, and the signal f(x) which is formed from the image. The dj (x) 
that is generated by this operation has the characteristic of rippling violently when 
the local periodicity of f(x) and the support of the mother wavelet \(/(x) (the interval 
at which \j/ is not 0) are roughly in coincidence (see FIG. 10). 

15 In the treatment shown in FIG. 9 and FIG. 10, the level j is changed in three 

stages for ease of explanation and understanding; however, theoretically, when the 
level j is changed in eight stages, the following takes place, 

2 8 " 1 = 2 7 = 128 -(17) 

and tree ring widths of over 100-fold can be addressed. In our experiment, it was 
20 verified that setting the level j in eight stages enabled us to accommodate tree ring 
widths of 0.1 mm to over 1 cm. Therefore, the number of stages of level j can be 
set at any number according to the characteristics of the object to be measured. 
9. Calculation of Squared Mean for Interval (Step SI 6) 

The squared mean for interval is calculated as shown in the following 
25 formula (18) for the dj (x) generated by the "8. Wavelet Conversion" (see FIG. 1 1). 
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gj (x) = I'-'a- 1 °\d j (ktfdk ... (18) 

Here, a is a constant that is determined by the support of the mother wavelet 

\|/(x). 

In the formula (18) the interval for which the squared mean is calculated is 
from x - 2~ } a to x + 2~ j a, but this is not the only possibility, and it can be defined as 
in the following formula (19). 



= 2 " , p;'rK( t )f* -d9) 



That is, since the pj is set in accordance with scale level j so that the pj 
becomes larger as the scale level j becomes lower, it is acceptable to make the 
10 calculation of the squared mean for interval for narrow intervals [x - pj, x + pj] for 
dj (x) of high frequency at high level j, and for broad intervals [x - pj, x + pj] for dj f 
of low frequency at low level j. The integration interval of x - 2' } a through to x + 
2~ j a, shown in the formula (18), is nothing but a standard integration interval. 
When the integration interval [x - pj, x + pj] is made narrower than this, it becomes 
15 more sensitive to local density changes, but more susceptible to the effects of noise. 
On the other hand, when the integration interval [x - pj, x + pj] is made broader than 
this, it loses its sensitivity to local density changes, but tends to become less 
susceptible to noise and more stable. Therefore, by adjusting the integration 
interval [x - pj 5 x + Pj] from the above standard state, it is possible to adjust the 
20 sensitivity of detection. 

10. Determination of the Dominant Level (Step SI 7) 
At a given point x, the level j at which thegj (x) defined in the formula (18) or 
formula (19) becomes maximum is the dominant level ja. Hence all points on the 
measurement line are made to have one dominant level (see FIG. 12). While the 
25 dominant level ja tends to become lower at portions where the tree ring width is 
wide, it tends to become higher at portions where the tree ring width is narrow or in 
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the vicinity of late wood where the signal changes markedly. 

1 1 . Determination of Maximum Value within the Interval (in the case of 
edge signal, the minimum value within the interval) (Steps SI 8 through S20) 
As shown in FIGS. 13A, 13B 5 and 13C, 
5 (a) When f(x) is a peak signal, when the value of f(x) is equivalent to the 

maximum value f ma x(x) of f(x) in the interval [x - 2" jd a, x + 2" jd a], then point x is 
determined as the point of maximum density within the tree ring layer. 

(b) When f(x) is an edge signal, when the value of f(x) is equivalent to the 
minimum value f m i n (x) of f(x) in the interval [x - 2" jd a, x + 2" jd a], then point x is 
10 determined as the late wood end within the tree ring layer. 

By conducting this operation throughout the entire range, it is possible to 
detect all tree rings on the measurement line. 

It should be noted that while in the above (a) and (b), the b in the 
abovementioned scaling parameter b* is 2, it can also be illustrated as follows. 
15 (c) When f(x) is a peak signal, when the value of f(x) is equivalent to the 

maximum value f ma x(x) of f(x) in the interval [x - b' jd , x + b~ jd ], then point x is 
determined as the point of maximum density within the tree ring layer. 

(d) When f(x) is an edge signal, when the value of f(x) is equivalent to the 
minimum value f m i n (x) of f(x) in the interval [x - b" jd , x + b" jd ], then point x is 

20 determined as the late wood end within the tree ring layer. 

To further generalize this, it can be illustrated as follows using the constant 
qjd that corresponds to the dominant level jd. 

(e) When f(x) is a peak signal, when the value of f(x) is equivalent to the 
maximum value f ma x(x) of f(x) in the interval [x - qjd, x + qjd], then point x is 

25 determined as the point of maximum density within the tree ring layer. 

(f) When f(x) is an edge signal, when the value of f(x) is equivalent to the 
minimum value f m i n (x) of f(x) in the interval [x - qjd, x + qjd], then point x is 
determined as the late wood end within the tree ring layer. 
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12. Cross-referencing of Main Measurement Line and Subordinate 
Measurement Lines (Step S21) 

When the tree ring on the main measurement line 1 established in the above 
"2. Designation of Measurement Site" is unclear, efforts to detect the maximum 
5 density point and late wood end of each layer according to the procedures in the 
above M 3. Acquisition of Pixel Information" through to "11. Determination of 
Maximum Value within the Interval (in the case of edge signal, the minimum value 
within the interval)" may not result in detection. In addition, partial unclarities of 
tree ring layers often occur randomly. This situation shall be referred to as non- 
10 detection. 

In addition, when there are points on the main measurement line 1 where the 
wood tissue is not uniform, or where there are density changes outside of the natural 
tree rings, such as cracks or polish marks (hereinafter referred to as noise), points 
other than the maximum density point or late wood end may be detected. This 

15 situation shall be referred to as erroneous detection. 

Normally, tree rings occur in concentric circles at the cross section and in 
parallel patterns at the radial section, while noise occurs in uncertain patterns. 
Taking advantage of these differences, verification is conducted as to whether a 
detection should be established as a tree ring, cross-referencing with the results of 

20 detection from the subordinate measurement lines configured in "2. Designation of 
Measurement Site" (the results of detection from each of the subordinate 
measurement lines 2 to 5 according to the procedure in the "3. Acquisition of Pixel 
Information" through to "11. Determination of Maximum Value within the Interval 
(in the case of edge signal, the minimum value within the interval)") in order to 

25 ensure the output of correct detection results even in the event of non-detection or 
erroneous detection on the main measurement line. The specific procedures for the 
above are described below. 

In the actual program, in order to conduct cross-referencing efficiently, the 
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positions and measurement line numbers of detection points on all main and 
subordinate measurement lines are recorded in memory as establishment candidate 
points. 

In the cross section, tree rings usually present a concentric circular pattern, 
5 but when considering an extremely small width, they can be regarded as being 
approximately parallel. The tree rings in the radial section are in a parallel pattern. 
However, since neither of them are in a complete parallel state, where they intersect 
orthogonally with the direction of the measurement line, there may be cases where 
there is a slight discrepancy between the detection position on the main measurement 

10 line and the detection position on the subordinate measurement lines. In order to 
absorb the error of the detection points, which should rightfully be the same tree ring 
layer, attributable to the different configuration position of the main or subordinate 
measurement lines, that is, in order to determine whether the detection point is within 
a roughly identical distance from the starting point of each measurement line, a width 

15 h, which can be regarded as being identical, is configured. 

Next, a definition is made as to how many detection points must occur out of 
the multiple main and subordinate measurement lines in order for that detection point 
to be established as a correct detection point. In our experiment, we obtained the 
best results with domestically produced Japanese cypress (Hinoki) cypress when 

20 detection points occurred in at least two out of the three measurement lines (the main 
measurement line and two subordinate measurement lines) in the cross section, and 
when they occurred in at least three out of the five measurement lines (the main 
measurement line and four subordinate measurement lines) in the radial section. 
More generally, the correct detection point can be defined as the detection point that 

25 has at least m detection points out of the n number of main and subordinate 
measurement lines within the width ±h that can be regarded as being identical. 

In accordance with the above rule, all the detection points on the main 
measurement line are subjected to analysis to establish those with detection points in 



at least m out of the n number of measurement lines within the width ±h as the 
correct tree ring (the maximum density point within each layer, or the late wood end 
within each layer) (see FIGS. 14A and 14B). Hence, 

(1) when the detection point on the main measurement line is the correct 
5 tree ring, then it is established as correct by cross-referencing. For instance, in the 

example in FIG. 15, from detection points ai through to a 3 there is a detection point 
on the main measurement line as well as detection points on at least three of the five 
measurement lines, and therefore the detection point is established. 

(2) when the detection point on the main measurement line is an erroneous 
10 detection, it is eliminated by cross-referencing. For instance, in the example in FIG. 

15, since the detection point a 5 has points in less than three locations out of the five 
measurement lines, it is not established but is subjected to the following "13. Cross- 
referencing between Subordinate Measurement Lines". 

(3) when there is a non-detection on the main measurement line (see 

15 detection point a 4 in FIG. 15), this situation is not considered in this process, and 
therefore the following "13. Cross-referencing between Subordinate Measurement 
Lines" is conducted. 

In "13. Cross-referencing between Subordinate Measurement Lines," as 
shown in FIG. 16, because detection point a 4 has detection points in at least three 

20 locations, it is established as priority 1 during the search. However, since detection 
point a 5 has points in less than three locations out of the five measurement lines, it is 
not established. 

In the stage shown in FIG. 17, search points aj through to a 4 have already 
been established, while on the other hand, detection point a 5 is not established as a 
25 detection point, and cross-referencing is finished. 

In order to conduct the following "13. Cross-referencing between 
Subordinate Measurement Lines" efficiently, those points that were established as 
correct tree rings in the process of "12. Cross-referencing of Main Measurement Line 
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and Subordinate Measurement Lines" are deleted in sequence from the establishment 
candidate points. 

13. Cross-referencing between Subordinate Measurement Lines (Step S22) 
The multiple subordinate measurement lines 2 to 5 are prioritized in advance. 

5 Usually, the closer it is to the main measurement line 1, the higher it is ranked in 
priority. In prioritizing, it must be ensured that no two measurement lines have the 
same rank. 

First, all the detection points remaining as establishment candidate points on 
the highest ranked subordinate measurement line are subjected to analysis to 

10 establish them as correct tree rings in accordance with the rule similar to that of "12. 
Cross-referencing of Main Measurement Line and Subordinate Measurement Lines" 
(see detection point a 4 in FIG. 16). Those points that were established as correct 
tree rings in this process are deleted in sequence from the establishment candidate 
points. A similar process is conducted in sequence from the highest priority down. 

15 This series of processes are finished when there are no more establishment candidate 
points, or when the process has been conducted on the lowest ranked subordinate 
measurement line (see FIG. 17). In reality, when the rule is "at least m out of n 
number of lines", it is meaningless to conduct analysis on the lowest ranked m - 1 
number of lines including the lowest ranked line; therefore, it is acceptable to cut off 

20 analysis at this rank. 

As described above, even in cases where there is a non-detection on the main 
measurement line in the process of "12. Cross-referencing of Main Measurement 
Line and Subordinate Measurement Lines", it is established as a correct tree ring by 
conducting "13. Cross-referencing between Subordinate Measurement Lines" (see 

25 detection point a 4 in FIG. 15 and FIG. 16). 

14. Determination of Tree Ring Location(Step S23) 

The established tree ring points of "12. Cross-referencing of Main 
Measurement Line and Subordinate Measurement Lines" and "13. Cross-referencing 



between Subordinate Measurement Lines" are consolidated to make the tree ring 
point (the maximum density point within each layer, or the late wood end within 
each layer) (see FIGS. 14A and 14B). The establishment points established on the 
subordinate measurement lines in the process of the above "13. Cross-referencing 
5 between Subordinate Measurement Lines" are projected onto the main measurement 
line at the time of consolidation. 

15. Measurement of Tree Ring Width and Output of Results (Steps S24 and 

S25) 

Because all the information on established tree ring points is amassed on the 
10 main measurement line 1, the number of pixels of each layer is counted from the 
coordinate value of each established tree ring point on the main measurement line. 
Because the size per pixel is determined by the image resolution set in the above "1. 
Acquisition of Tree Ring Image", 

(Tree ring width) = (^ r m , b a ye? f ^'"l X < size ? er pixel) 

15 allows for the calculation of the tree ring width of each layer. By output of this 

result, the entire process is completed. 

Incidentally, tree ring width is not the only data that is useful in 

dendrochronology. Other important data include maximum density within the layer, 

minimum density within the layer, and early wood/late wood ratio. 
20 The effect of the method for the measurement of waveform signal feature 

points by the process described above is shown in Table 1 and Table 2. 
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Results of detection of maximum density point using dens i ty s i gna 1 
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Table 2 

Results of detection of late wood end using differential signal 
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Table 1 shows the results of detection of maximum density points when the above 
mentioned density point is used as a waveform signal, compared against detection 
5 results using the conventional method. Table 2 shows the results of detection of 
maximum late wood end when the above mentioned differential signal is used as a 
waveform signal, compared against detection results using the conventional method. 
In Table 1 and Table 2, the new method according to the present invention is 
separated into wavelet conversion only (cross-referencing of multiple measurement 
10 lines omitted), cross-referencing of multiple, measurement lines only (no wavelet 
conversion), and the combined use of wavelet conversion and cross-referencing of 
multiple measurement lines. 



Table 3 

Definitions: Rates of detection, erroneous detection and 
non-detection 





Detect i on 


Detected 


Not detected 


Input 


Tree ring 


Ss 


S N 


Not tree r ing 


Ns 


Nn 



30 



The rate of detection, rate of erroneous detection, and rate of non-detection 
shown in the above Table 1 and Table 2 are defined by the method shown in the 
following formulas (21), (22), and (23) from the relationship of Ss, Sn, Ns, and Nn 
shown in the above Table 3. 

Rate of detection = Ss/ (Ss+Sn) •••(21) 
Rate of erroneous detection = Ns/ (Ss+Ns) -"(22) 
Rate of non-detection = Sn/ (Ss+Sn) •••(23) 

5 

As can be understood from the above Table 1 and Table 2, it was verified 
that in all cases, using the measurement method of the present invention resulted in 
increased rate of detection and decreased rate of erroneous detection and rate of non- 
detection, and it was proven that detection performance improved overall. 

10 The above has been a detailed description of measuring tree ring width, and 

if the process up to and including "14. Determination of Tree Ring Location" can be 
conducted with accuracy, it is easily possible to calculate, from peak signals and 
edge signals, characteristic quantities other than these tree ring widths. 

In addition, the measurement method and measurement equipment according 

15 to the present invention are not limited in their use to tree ring measurement of wood 
specimens, but can also be applied to such uses as measurement of fingerprints, voice 
prints, and retina patterns for authentication purposes, measurement of wiring 
patterns in electronic components, and measurement of such biological signals as 
brain waves. 

20 INDUSTRIAL APPLICABILITY 

As is clear from the above description, the method and equipment for the 
measurement of waveform signal feature points according to the present invention 
allows for accurate detection of the feature points of peak waveform signal even 
when the waveform signal is irregular in its distance between feature points; 

25 therefore, it is possible to markedly improve the detection performance of tree rings 
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etc. in wood specimens. 

That is, before incorporating wavelet conversion into the measurement 
method, tree ring detection was conducted by simple binary coding wherein a 
threshold value was established on f(x) and a given value was distinguished by being 
5 either over or under the threshold value, therefore having the drawbacks of being 
extremely vulnerable to noise and/or indistinct tree rings, and of the threshold value 
tending to be dependent on differences between individual specimens; however, by 
incorporating the measurement method using the wavelet conversion according to the 
present invention, it is possible to conduct tree ring detection while automatically 

10 adapting to the fineness (or coarseness) of the tree rings, or the clarity (or 
indistinctness) of the tree rings, as each situation arises, which contributes to the 
improvement in detection performance. 

In addition, since this detection performance does not require the setting of a 
threshold value, it has a robust aspect that is not easily swayed by differences 

15 between individual specimens, and in this respect also is extremely innovative and 
effective. 



